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We show how the phases of interacting particles in topological flat bands, known as fractional 
Chern insulators, can be adiabatically connected to incompressible fractional quantum Hall liquids in 
the lowest Landau- level of an externally applied magnetic field. Unlike previous evidence suggesting 
the similarity of these systems, our approach enables a formal proof of the equality of their topological 
orders, and furthermore this proof robustly extends to the thermodynamic limit. We achieve this 
result using the hybrid Wannier orbital basis proposed by Qi [Phys. Rev. Lett. 107, 126803 (2011)] 
in order to construct interpolation Hamiltonians that provide continuous deformations between the 
two models. We illustrate the validity of our approach for the groundstate of bosons in the half 
filled Chern band of the Haldane model, showing that it is adiabatically connected to the v — 1/2 
Laughlin state of bosons in the continuum fractional quantum Hall problem. 



Owing to the recent discovery of topological insulators 
PQ, there is now hope to realize materials which mani- 
fest Haldane' s vision of a quantum Hall effect without 
external magnetic fields [2j [3]. Several proposals have 
extended this concept to fractional quantum Hall (FQH) 
liquids that could be realized in topologically non-trivial 
bands which are also flat |3H5]. A similar mechanism 
was proposed to simulate the effect of strong magnetic 
fields in cold atomic gases [9j [10] . Flat bands with non- 
zero Chern number [TTJ [12] provide an avenue to realize 
strongly correlated states at high temperatures [4], rais- 
ing prospects for stabilizing exotic non-abelian phases re- 
quired to build topological quantum computers [15] . 

Numerical works seeking evidence for incompressible 
quantum liquids in topological flat bands have focused 
on spin polarized models breaking time-reversal symme- 
try [T4HT6] , which were baptized as fractional Chern in- 
sulators (FCI) [16]. Signatures for the topological na- 
ture of their ground states include their spectral flow and 
groundstate degeneracies p~4j [15] and the analysis of the 
entanglement spectra. The latter reveal a counting of 
excitations matching that of FQH states at the corre- 
sponding band filling, e.g. the Laughlin and Moore-Read 
[T6] or Jain states [17] . One drawback is that these data 
can be acquired only for finite size systems. 

Current insights in the analytic theory of fractional 
Chern insulators rely on the analysis of the projected 
density operator algebra p~8j [19] or of emergent sym- 
metries in the exact many-body spectrum [20] . There 
are several proposals for constructing FCI wave functions 
[2TH26] , and some overlaps were calculated for fermionic 
systems [26]. However, the understanding of the many- 
body ground states of FCIs cannot yet pride itself with 
an achievement similar to the celebrated accuracy of an- 
alytical wave functions for FQH states [27H29] . 

In this manuscript, we provide a formal proof that FCIs 



are in the same universality class as FQH states. We 
base our argument on Qi's proposal [21] of a mapping 
between FQH and FCI wavefunctions. By representing 
both FCI and FQH Hamiltonians in a Hilbert space with 
the same structure, we are able to study a class of su- 
perposition Hamiltonians that extrapolate smoothly be- 
tween these systems. Taking advantage of this construc- 
tion, we demonstrate that the many-body ground states 
of bosons in a half filled LLL and the topological fiat 
band of the Haldane model are adiabatically connected, 
proving formally that these phases have the same type of 
topological order in the thermodynamic limit. Our strat- 
egy can be employed generally to identify incompressible 
quantum liquids in Chern bands, including bands with 
Chern number \C\ > 1. 

Let us briefly comment on the case of lattice FQH 
states [30H33] , which are simultaneously FCIs in the 
Chern bands of the Hofstadter butterfly [34] given that 
flux can be gauged away. Nonetheless, these systems can 
be taken to the limit of the continuum FQHE [35 . Hof- 
stadter bands can be realized in cold atomic gases [36] , 
which may also provide the most promising avenue for 
realising FCIs since topological flat bands require fine- 
tuned parameters that are efficiently controlled in these 
systems HUJ [36] ES] • 

We first establish our notations for the description 
of fractional Chern insulators. We consider finite two- 
dimensional lattices of iV ce n = L\ x L 2 unit cells, spanned 
by lattice vectors forming an opening angle 7, and we 
choose vi = sin(7)e^ + cos(7)e y and v 2 = e y . Lat- 
tice sites are located on sublattices a within the unit 
cell. In a finite system with periodic boundary con- 
ditions \I/(r + L{Vi) = e 1( ^\I/(r), the reciprocal lattice 
consists of discrete points k = J^fe + 2^)*^' wnere 
Gi = 2iie x /Li sin(7) and G 2 = 2tt[— cot^e^ + e y ]/L2 
and we consider a rhomboid fundamental region with 



2 



qi = 0, ...,L i - 1. 

For the moment, we consider the Hamiltonian of the 
infinite system in its momentum space representation 
H = Ek4, a ^(k)cM, where ft(k) yields a topologi- 
cally non-trivial flat band. Let the Bloch functions u 
and eigenenergies e be determined by the correspond- 
ing eigenvalue equation h a p(\t)u 7 p = e n (k)i£™, intro- 
ducing the band index n, and with the normalization 
|ii™(k)| 2 = 1. In the following, we will denote the 
eigenstates as |n, k) = u^k)^ a |vac). The ensuing 
Berry connection *4(n, k) = — i ^ a (k)V ^u^k) is a 
gauge dependent quantity. Physical observables such as 
the Berry curvature B(k) = Vk A A.{k) and the resulting 
Chern index C = ^ j BZ <i 2 kS(k) are gauge invariant. 

We now review Qi's proposal for mapping FCIs onto 
FQH states by a construction of Wannier states within 
a topological flat band [21] . In his approach, we formu- 
late hybrid Wannier wavefunctions \W(x,k y )) that are 
localized only along the x-axis, while retaining trans- 
lational invariance with a well defined momentum pro- 
jection k y onto the ^/-direction [39]. We can think of 
these states as the simultaneous eigenstates of momen- 
tum P y and the band-projected position operator X cg = 



\im n 



1 9 n 



, satisfying 



X c ° \W(x, k y )) = [x- 0(ky)/2n] \W(x, k y )) . (1) 



We adopt the explicit construction of the Wannier 
states in terms of the momentum eigenstates of h(k) 
presented in Refs. [2TJ [40], given by \W(x,k y )) = 
Y.ll=ofk x ky \ n = d k x,k y )), for x = 0, - 1, with 



(2) 



This expression is related to a simple Fourier trans- 
form of the momentum eigenstates by additionally tak- 
ing into account parallel transport of the phase along 
k x according to the Berry connection A x . The polar- 
ization 0(k y ) = Jq^ A x (p Xj k y )dp x is required to ensure 

periodicity of the state in k x , enforcing f^ y = /^+ 2 7r- 
The relative phases x(ky) °f the Wannier states repre- 
sent a gauge freedom of the theory, while the relative 
phase of Bloch functions at the same k y is fixed by the 
Berry connection in (J2|. We take the particular choice 

X(k y ) = exp[-if£ y Ay(0,py)dp y + i^ J 27r Ay(0,py)dp y ], 
as suggested in [40] . 

In finite-size systems, we adapt the construction ([2| 
straightforwardly. Given the Bloch functions u^(k) 
on reciprocal lattice points in the fundamental region, 
we choose a gauge that is consistent with the peri- 
odicity of momentum space, i.e. ^(k -f Z^Gi) = 
i^(k). A discretized version of the Berry connec- 
tion of band n is then computed as A™(qi,q2) = 
3 log [<*(<?!, g 2 )<(<2i + 1,22)]. The integral of the 
Berry connection is discretized as f x A x (p x ,k y )dp x — >• 




FIG. 1. (a) Geometry and hopping terms in the Haldane- 
Model: the fundamental unit cell has two inequivalent sites 
A and B. Second nearest neighbor interactions are com- 
plex, with arrows indicating the direction of a positive hop- 
ping phase (j) rr f . (b) Berry curvature for the Haldane model, 
(c) Corresponding expectation value of the position operator 
(X cg ) for Wannier states according to Eq. (M). 



j^gife) A™(qi,q2), and mutatis mutandis for Ay(qi,q2). 
This resolution of ([2| yields a unitary transformation of 
the original single particle basis. As angles, the values of 
A x (qi,q2) are defined modulo 27r, and we ensure that the 
shift in x position Ax = 0(k y ) /2tt satisfies < Ax < 1. 
The Wannier states can thus be brought into an order of 
increasing centre of mass position (X cg ) by using a single 
linearized momentum index J [21 relating to the param- 
eters of the Wannier state by K y — k y + 2ttx = 271-J/L2, 
with J = 0, . . . , iVceii — 1. 

To describe the fractional quantum Hall problem of 
particles in the lowest Landau-level of a magnetic field, 
we adopt the Landau gauge A = —xBe yi such that 
our oblique simulation cell is pierced by = N ce \\ 
flux quanta, with L1V1 x L2V2 = 2tt£qN ( / ) . The peri- 



odic Landau- level orbitals (f)j(x,y), j = 0, ...,iV^ — 1 
[4TJ [42] are chosen with definite momenta k y = 27T//Z/2, 
and achieve their maximum amplitude at (x) = kyi^. 

For the remainder of this manuscript, we choose a spe- 
cific flat band model to perform a quantitative assessment 
of the Wannier representation: we work with the Haldane 
model [2], defined on the lattice shown in Fig. [TJa), and 
choose parameters yielding a nearly flat C — 1 band: 
t x = 1, t 2 = 0.60, t 3 = -0.58 and <j> = 0.4tt [15 . The 
resulting Berry curvature, shown in Fig.JIJb), is inhomo- 
geneous. The definition implies that 



_d_ 

dk v 



(X C9 )U 



1 d0(k y ) 
2tt dk y 



p2tt 

L • 



B(p x ,k y )dp x , (3) 



i.e., the k y dependency of the integrated Berry curvature 
translates into a non-linear evolution of the centre of mass 
position for the Wannier states, shown for the Haldane 
model in Fig. [TJc). By contrast, the lowest Landau-level 
has a constant Berry curvature, yielding linear behaviour. 

Having defined a single particle basis {<j)j} character- 
ized by a single linear (linearized) momentum index j 
(J) for the lowest Landau-level (FCI), respectively, we 
can compare the structure of their interaction Hamil- 
tonians by evaluating matrix elements. Formally, two- 
body interactions can be written in the generic form H = 



£{i<} v hnhn^h^ c ^ with matrix elements v hnhn 
given by the projection to the lowest band p~6j|41]. 

We focus on contact interactions for bosons, as this 
case has a straighforward interpretation for the contin- 
uum and on lattices. To treat the FCI case, we flat- 
ten the residual dispersion of the topological band (en- 
suring that the Wannier states are energy eigenstates) . 
Upon comparison, we find that the matrix elements for 
the FCI Hamiltonian differ from the FQH case in two 
aspects: The first issue concerns momentum conserva- 
tion. For the FQHE, the momentum of the Landau- 
gauge k y is conserved in scattering processes, and hence 
yf 1 < ]fj 3 j 4 oc £j 1+ j 2j3+ j 4 . In a topological flat band, mo- 
menta qi and q 2 are conserved separately while the lin- 
earized momentum index of the Wannier states J is con- 
served only modulo L 2 , i.e., Vj^j A ex &™j 2 j\ + j A . In 
Fig.|2|a), we illustrate the magnitude of matrix elements 
for a small system. The figure clearly shows the block- 
diagonal structure for the FQH, reflecting full momentum 
conservation, while the FCI Hamiltonian has several off- 
block-diagonal entries. Nevertheless, both matrices are 
similar in that the entries of largest magnitude are lo- 
cated at the same positions. The second difference lies 
in the translational invariance of the matrix elements. 
For the FQHE, the amplitude for scattering processes is 
invariant under translations in momentum space (or ef- 
fectively in real-space, given that (x) oc k y ). For the 
FCI on the other hand, the non-linear depencency of 
(X C9 (K y )} imprints a variation of the matrix elements 
with periodicity L 2 . This effect is illustrated for sev- 
eral nearest neighbor interactions in Fig. [2|b,c). Given 
these two qualitative differences - momentum conserva- 
tion and translational invariance of the matrix elements 
- the FCI Hamiltonian in the Wannier basis cannot have 
eigenstates that are identical to those the corresponding 
FQH problem, as had been conjectured in Ref. [21 . 

As a next step, we evaluate the similarity of the wave- 
functions for the FCI and FQH problems in terms of their 
overlap when written in the Wannier and Landau-gauge 
bases, respectively. We analyze the case of a half filled 
band, or v = 1/2, for systems with TV = 6, 8, and 10 
bosons, on lattices of several aspect ratios. For the cor- 
responding FQH problem, we choose a simulation cell 
with the same geometric features, namely a torus with 
an aspect ratio given by R = L 2 /L\ and opening angle 
7 = 7r/3 to match the hexagonal lattice underlying the 
Haldane model. The FCI Hamiltonian in the Wannier ba- 
sis is diagonalized in the Fock spaces for total linearized 
momenta J to t = ^2n=i ^ n m °d L 2 . The Hilbert space for 
the FQHE has full translational symmetry and segments 
into blocks with total momentum j^ ot = 0, . . . , — 1, 
(where = N ce \\). Accordingly, each FCI eigenstate in 
sector J to t can have overlap with several sectors j^ ot of 
the FQH problem satisfying [j^ ot mod L 2 ] = Jtot- In ad- 
dition, the Laughlin state [27 which is the exact ground- 
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FIG. 2. (a) Magnitudes of matrix elements for the delta 
interaction between particles between pairs of incoming (out- 
going) momenta — [( Ji + J2) mod iVceii] in a finite size 
system of N — 6 particles for the FQHE on the torus with 
— 12 (left) and for the FCI in the lowest band of the Hal- 
dane model for a 3 x 4 lattice (right), (b) Schematic showing 
some short range interactions, including a 'squeezing' process 
Vn as well as two diagonal interaction terms V20, Vbo- Panel 
(c) shows how the magnitude of these processes depends on 
the centre of mass position for the FCI (solid) as compared 
to FQH (dashed). 



state of contact interactions in the LLL at filling v = 1/2 
has a twofold topological degeneracy cZqs = 2. Hence, we 
calculate the total groundstate overlap O as an average 
for both groundstates |^^ CI ), taking into account pro- 
jections VjT onto sectors with torus groundstates 

yielding O = ^ £f G s Li K^t^I*^ 1 )! 2 . For our 
largest system, N = 10 particles on a L\ x L2 = 4 x 5 
lattice with a Hilbert space of d = 4 x 10 6 , we find a 
clear gap above two low lying groundstates that yield a 
total overlap of O = 0.826. This value correspond to 
Qi's gauge-choice of the Wannier states [40 . In addition, 
we have run numerical optimizations of the phases x(fc 2/ ) 
maximizing (9, and have found minor changes < 1% in 
the overall result. Thus, we report overlaps conforming 
with the initial gauge choice [43]. The FCI states have a 
total weight W = 0.892 within the momentum sectors of 
torus groundstates, establishing an upper bound for the 
overlap. The 'leakage' of weight outside the groundstate 
sectors results from off-block diagonal entries in the FCI 
Hamiltonian, and is independent of the gauge choice. 

For the Haldane model, we conclude that the Wannier 
construction yields non-trivial overlaps with the eigen- 
states on the torus. However, in light of our results, Qi's 
construction does not yield satisfactorily accurate trial 
wavefunctions. Nevertheless, the Wannier construction 
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FIG. 3. (a) Overlap of the groundstate manifolds of 7-L(k) and 
H FQH (see text) for bosons at v = 1/2. (b) Total weight in 
torus groundstate subspace. (c) Spectrum for a system with 
N = 10 particles, along a path adiabatically connecting a 
continuum problem on the torus to the FCI Haldane model 
on a lattice of 4 x 5 unit cells. The insets show the momentum- 
resolved spectrum for the torus (left) and the pure FCI system 
(right), (d) Gap for several systems of different sizes and 
aspect ratios. Inset: finite size scaling of the gap. 



allows us to formulate the FCI and FQH problems in 
a Hilbert space of the same structure, making it conve- 
nient not only to calculate overlaps but also to construct 
an adiabatic continuation between them. Hence, we can 
formulate a superposition of both Hamiltonians and an- 
alyze its spectrum at any intermediate value of an inter- 
polation parameter n : 

H(k) = kVtfbU fg1 (U) + (1 - K )H FQH (Vo = 1), (4) 

where Vtfb denotes the flattening of the topological flat 
band. In order to fix the relative energy scales in the 
two problems, we analyze the magnitude of the gap and 
choose a value of U that equalizes its numerical value for 
n = and k = 1, respectively. As we find little scaling 
of the gap with the system size (see below), we choose a 
single value of U = 0.2649 throughout our study. With 
the definition Q of the adiabatic continuation, we can 
analyze the overlap and leakage of the groundstate wave- 
functions as a function of n. Fig. [3ja) summarizes our 
results for several system sizes, showing how the overlap 
drops with increasing system size. 

We now examine whether the Laughlin state on the 
torus is adiabatically connected to the groundstate of the 
Haldane model with delta interactions. The properties 
of topologically ordered phases are conserved under the 
variation of system parameters as long as the ground- 
state manifold is protected by a finite gap A. Hence, we 
evaluate the spectrum of the class of Hamiltonians Q 
as a function of as displayed in Fig. [3^c) for N = 10 
particles. The spectrum clearly shows a twofold degen- 
erate groundstate, which is well separated by a gap from 
a continuum of excited states at higher energy. To sur- 
vey finite size scaling, we report A for different lattice 
geometries in Fig. [3^d). The magnitude of A is found 
to be weakly dependent on the interpolation parameter 
n. Furthermore, it also depends weakly on system size. 
The inset shows the finite size scaling of A for different 
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FIG. 4. Entanglement spectrum (a) and entanglement gap 
(b) for Na = N/2, in analogy of the energy spectrum shown 
in Figged). 



adiabatic continuation parameters, clearly revealing that 
the gap remains open in the thermodynamic limit. Hence 
we confirm that the bosonic Laughlin state at v = 1/2 is 
adiabatically connected to the ground state of the Hal- 
dane model, which firmly establishes that they are in the 
same universality class. 

Finally, given its importance for identifying incom- 
pressible states in FCI models [I6j [44J , we consider the 
entanglement spectrum of the ground states along the 
trajectory < k < 1. We evaluate the particle entan- 
glement spectrum that encodes the number of quasihole 
excitations above the groundstate [45]. We find that the 
count of entanglement eigenstates below the entangle- 
ment gap A^ remains conserved at all values of n. Our 
data shows the occurrence of isolated, low-lying entangle- 
ment eigenstates in proximity of the universal eigenstates 
at low Schmidt eigenvalues (see Fig.|4|, however, we have 
not observed any closure of the gap. The magnitudes of 
the gap above the universal states are summarized in 
Fig. |4^b), confirming that the gap remains open, but the 
data do not allow a concise analysis of finite size effects. 
Nonetheless, our data are consistent with an extended 
adiabatic continuity in terms of the entanglement gap. 

In conclusion, we have established an approach to show 
that fractional Chern insulators are adiabatically con- 
nected to fractional quantum Hall states. Our technique 
uses Qi's construction of hybrid Wannier orbitals, and 
extends to the thermodynamic limit by a robust extrap- 
olation procedure. Specifically, we have used this con- 
cept to prove that the FCI ground state of bosons in the 
half filled Chern band of the Haldane model is in the 
same universality class as the Laughlin wavefunction at 
v = l/2. 
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edges support from the Leverhulme Trust under grant 
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